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ABSTRACT 

We report the serendipitous detection of a very bright, very nearby microlensing event. In late 
October 2006, an otherwise unremarkable AO star at a distance ~ 1 kpc (GSC 3656-1328) brightened 
achromatically by a factor of nearly 40 over the span of several days and then decayed in an apparently 
symmetrical way. We present a light curve of the event based on optical photometry from the Center 
for Backyard Astrophysics and the All Sky Automated Survey, as well as near-infrared photometry 
from the Peters Automated Infrared Imaging Telescope. This light curve is well-fit by a generic 
microlensing model. We also report optical spectra and Swift X-ray and UV observations that are 
consistent with the microlensing interpretation. We discuss and reject alternative explanations for this 
variability. The lens star is probably a low-mass star or brown dwarf, with a relatively high proper 
motion of > 20 mas yr , and may be visible using precise optical/infrared imaging taken several 
years from now. A modest, all-sky survey telescope could detect ~ 10 such events per year, which 
would enable searches for very low-mass planetary companions to relatively nearby stars. 
Subject headings: stars: individual (GSC 3656-1328); gravitational lensing 



1. INTRODUCTION 



At the very moment in 1936 that he introduced 17 the 
concept of gravitational microlensing of one star by an- 
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other closely aligned along its line of sight, Einstein fa- 
mously dismissed its practical significance. Noting that 
the characteristic scale (what we now call the "Einstein 
radius") was extremely small, he concluded that "there 
is no great chance of observing this phenomenon, even 
if dazzling by the light of the much nearer star is dis- 
regarded." (|Einsteinlll936D . It is easily shown that the 
optical depth to microlensing (the probability that any 
given star lies projected within the Einstein radius of an- 
other) is only about r ~ 10~ 8 among the V < 12 stars 
that were typically cataloged in Einstein's day. As there 
are only a few million such stars altogether, the probabil- 
ity that any of these are microlensed at any given time is 
much less than one (|Nemiron1ll998[ ). Such a calculation 
may have influenced Einstein to resist the determined ef- 
forts by Hungarian engineer R. W. Mandl to get Einstein 
to publish his microlensing formulae and perhaps also to 
compose a letter to the editor of Science to "thank you 
for your cooperation with the little publication, which 
Mister Mandl squeezed out of me. It is of little va lue, 
but it makes the poor guy happy" (|Renn et al.lll997t ). 

Since Einstein's 1936 article, several authors h ave at- 
tempted to resurrec t the idea of microlensing (e.g. JLiebesI 
1964; Rcfsdal 1964). However, microlensing experiments 
did not actually get underway until the early 1990s 
(lAlcock et all 119931 : lAubourg et "all 119931 : lUdalski et "all 
1993). These experiments wer e motivated both by the 
suggestion of iPaczvhskil (|1986h to simultaneously mon- 
itor millions of stars in the dense star fields of nearby 
galaxies and by contemporaneous advances in technol- 
ogy that made such experiments feasible. 

To date, several thousand microlensing events have 
been discovered toward several lines of sight, including 

for the deflection of light by a massive body over 100 years before 
Einstein, although he did not c onsider the associated m agnification 
of the source by the lens. See Schneider et al. ( 1992) for a more 
thorough discussion of the history of gravitational lensing. 
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the Large and Small Magellenic Clouds (lAlcock et al 
1997[ iPalanque-Delabrouille et alJ 119981: 
2001 

200: 



as well as M31 (iPaulin-Henriksson et al 
Ide Jong et all 12004 lUglesich et~aTl 12004 
Calchi Novati et al r i2005f ). However, the vast ma- 



jority of ev ents have been detec t ed toward the Galac- 
tic bulge (lUdalski et all 120001: iThomas et ail 120051 : 
lHamadache et al.l I2006T) or fie lds in the Galacti c plane 
relatively close to the bulge (jDerue et al.1 1200 ID . The 
source stars of these events have generally been relatively 
faint, V > 15. In these cases, the optical depth is of 
order 10 -6 , i.e. 100 times higher than for the local stars 
that Einstein would have considered because the targets 
are about 10 times farther away. The larger distance 
makes the area of the Einstein ring about 10 times 
bigger and increases the column density of potential 
lenses by another factor of 10. 

In the intervening years since the first microlensing 
events were discovered, a few authors have revisited 
the idea o f detecting bright an d/or nearby microlens- 
ing events. ICollev fc Go tt ( 1995) argued that microlens- 
ing events visible to the naked eye are exceedingly rare, 
occurring at a rate of one per 40,000 yr for lensing by 
known stars. Microlensing of fainter stars is obviously 
more common; the ev ent rate for stars with V < 15 is 
roughly one per year (Nemirofl 1998; Han 2007). This 
calcul ation led iNemirofj (fl 998( 1 and iNemiroff fc RafertJ 
(1999) to point out that only small-aperture instru- 
ments are required to discover these brightest microlens- 
ing events, and that monitoring the entire sky down to 
V ~ 15 on relativel y shor t timescales would soon be 
feasible. Di Stefano (2005) considered in detail a re- 
lated idea of "mesolensing" : microlensing of background 
stars by nearly lenses with large angular Einstein rings 
and large proper motions. Mesolensing can be used to 
study the properties of nearby stars, and in particu- 
lar their planetary comp anions dDi Stefandl2?X)7ll200ST : 
iDi Stefano fc Night! 120081 . See SJfor further discussion 
of the potential of microlensing to discover and study 
nearby planetary systems. 

Although, given the typical source distances of ~ 
8 — 50 kpc for microlensing searches toward the Galac- 
tic bulge and the Magellanic Clouds, the majority of 
the microlenses giving rise to observed events have been 
located at distances of a kiloparscc or more, there is 
neverthele ss a low-probabili t y tail of more n earby mi- 
crolenses jGouidiLan GJH ETStefen^IlOOl). Indeed, 
there is one event for which the lens has been robustly lo- 
cated to within a ki l oparsec of the Sun , MACHO-LMC-5 
(| Alcock et all 120011 : iGoiild et al.ll2004 l. Another event, 
EROS2-LMC-8, is also likely due to a nearby microlens, 
although this has yet to be confirmed with follow-up ob- 
servatio ns (J. P. Beaulieu 2008, priv ate communication; 
see a lso iTisserand fc Milsztainll200"4l and iTisserand et al.l 
12007( 1. Thus, nearby lenses can be uncovered and studied 
in ongoing surveys. However, these surveys are not very 
efficient at discovering nearby microlenses; furthermore, 
the source stars of these events are necessarily faint, mak- 
ing detailed monitoring and follow-up observations diffi- 
cult. 

Here we report on a microlensing event of the type Ein- 
stein believed would never be observed: a magnification 
A = 40 event of the V ~ 11.4 A-type star GSC 3656- 
1328, located about 1 kpc from the Sun in the disk of 



the Milky Way. Although such events are indeed quite 
rare, microlensing of somewhat fainter nearby stars oc- 
curs with reasonable frequency. While this event was 
found serendipitously, we argue that with recent techno- 
logical advances it is now feasible to monitor the sky to 
deliberately and routinely detect these "domestic" mi- 
crolensing events. We propose a telescope design to ac- 
complish this, and we argue that it is possible to build 
two copies of such a telescope that could monitor the 
majority of the sky down to V ~ 16 at modest cost. 
Monitoring of the discovered microlensing events would 
enable a novel method to detect nearby planets, allow a 
search for dark objects in the Milky Way disk, and permit 
several days advance warning for potentially hazardous 
a steroid impacts. 

iFukui et alJ (|2007l ) also report on observations of the 
transient event in GSC 3656-1328. A subset of the data 
presented here is in common with theirs; however, the 
analyses were done completely independently. They also 
conclude that the brightening seen in GSC 3656-1328 is 
likely due to microlensing. 

2. OBSERVATIONS 

On 2006 October 31, A. Tago announced a sud- 
den brightening in the close vicinity of GS C 3656-1328 
(jNakano fc Tagoll2006t iNakano et al J 12006( 1. a V = 11.4 
A0V-A1V star 18 ~ 8° from the Galactic plane in Cas- 
siopeia with a distance of ~ 1 kpc. Confirmation fol- 
lowed a few hours later in the variable-star newsgroups; 
this established that the star was in fact the source of 
the brightening, and roughly measured the amount and 
timescale of the brightening (to V=7.5, rising in about a 
week) . 

We began our photometry with the small tele- 
scopes of the Center for B ackyard Astrophysics (CBA: 
ISkillman fc Pattersonl[l993( l on November 1 (=JD'=JD- 
2450000=4040.), and found the star at V = 8.9, falling 
rapidly and smoothly, with no additional variability at 
the few percent level. In the next few days, we ob- 
tained time-series and multicolor photometry with CBA 
telescopes, spectroscopy with the MDM 2.4 m tele- 
scope, aji^_a_p_clirted_5000 s X-ray ob servation with 
Swift (|Gehrels et al.ll2004 ) using the XRT (jBurrows et all 

2005) ; we also searched for X-ray outbursts over the 10-yr 
lifetime of the All-Sky Monitor aboard RXTE. The ba- 
sic conclusion from all of these observations was simple: 
the star had brightened by 4 mag, but with no change in 
spectrum or color (B—V = 0.2), no flickering, and no dis- 
cernible X-ray signature (with a Swift upper limit in the 
0.5-10 keV band of < 4x 10~ 13 ergs cm~ 2 s _1 ). A search 
of 400 photographic pl ates during 1964-1994 s howed no 
variability in the star (|Samus fc Antipinir2006h . and we 
found no X-ray outbursts from this position. 

This seemed unlike any known type of intrinsic vari- 
able star, and several reports suggested instead that 
GSC 3656-13 28 may have microlensed b y an intervening 
passi ng star (jMikolajewski et all l2006al fbl: iSpiegel et al.1 

2006) . After 2 weeks of data, our preliminary fit of 
the light curve to a microlens model (see $3]) seemed 
promising. However, there were only a few very uncer- 
tain measurements (with typical errors of ±0.4 mag) be- 

18 We summarize the properties of GSC 3656-1328 in Appendix 
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fore the peak, and identifying microlensing events based 
on falling-only light cur ves can be problematic (jSmithl 
120031 1 Afonso 'et al.|[200l) . Given the low a priori likeli- 
hood for such a nearby star to be microlensed, we eagerly 
sought additional data to confirm (or refute) the mi- 
crolensing interpretation. In particular, we sought data 
before the apparent peak on October 31, as well as data 
covering a wide baseline in wavelength, in order to as- 
certain whether the variation was both symmetric about 
the peak and achromatic, as would generally be expected 
for a short microlensing event of an isolated star. 

We were fortunate to find 15 images of the field in the V 
and I bands in the test runs of the northern (Hawaii) sta- 
tion of the All-Sky Automa ted Survey (AS A S; the south- 
ern station is described bv lPoim ahski 20041). These im- 
ages have errors of ±0.04 mag, cover (athough sparsely) 
the rise of the event, and fortuitously include one I- 
band measurement obtained almost exactly at the peak, 
judged by a microlensing fit to the remainder of the data 
set, as described below. Since ASAS-North had just be- 
gun operation, this was mighty lucky. 

As for wavelength coverage, we supplemented the CBA 
monitoring with four epochs of UV observations (on 
JD' = 4042, 4055, 4057, and 4067 ) using the UVOT cam- 
era on Swift (|Roming et al.|[2005h , and a long program of 
infrared monitoring from the Pet ers Automated Infr ared 
Imaging Telescope (PAIPJTEL, iBloom et all 120061 ) on 
Mount Hopkins. PAIPJTEL is a 1.3m telesc ope eq uipped 
with the former 2MASS (jSkrutskie et alJ [2006) cam- 
era that simultaneously images the near infrared JHK S 
bands. We obtained ~ 1500 7.8 s images in each band 
nearly every night from program start ( JD' = 4044) until 
the event returned to baseline. Since the variations were 
smooth, we binned the images into 51 separate epochs, 
and the magnitude of the source at each epoch was de- 
termined using differential aperture photometry against 
a set of reference stars from the 2MASS catalog. The 
Swift UV data from this bright source were strongly af- 
fected by coincidence losses, and to minimize these we 
used only the M2 band (centered on 2400A) and esti- 
mated the flux of the source by summing in an annulus 
centered on the source, with an inner radius of 12" and 
an outer radius of 20". 

We also obtained seven optical spectra using the CCDS 
instrument on the MDM 2.4m telescope on Kitt Peak, 
two near JD' ~ 4043.7, three near ~ 4046.7, and two 
near ~ 4057.7. The spectra cover 4000 - 6800A, with a 
resolution of ~ 15A FWHM, as measured from arc-lamp 
lines. They reveal an unremarkable early A star, showing 
prominent resolved (FWHM~ 25A) Balmer lines. We 
measured the equivalent widths (EWs) of Ha, H/3, and 
H7 in each spectrum by integrating over a window of 
±60A, fitting a third-order polynomial to estimate the 
local continuum. We estimated the uncertainties in the 
EW measurements by generating a series of mock spectra 
with identical noise properties, and measuring the EWs 
in these spectra in the same manner as the actual data. 
See Figure [TJ 

3. MICROLENSING INTERPRETATION 

Figure [2] shows the collected UV, optical, and near- 
IR photometry of GSC 3656-1328 during the transient 
variability. For all of these data except those from the 
CBA, the photometric uncertainties are due to the pho- 
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Fig. 1. — EW of the Ha (red circles), H/3 (blue triangles), 
and H7 (green circles) lines of GSC 3656-1328 as a function of 
JD-2450000. The measurements for H/3 have been offset by 2Afor 
clarity, as indicated. The broadband optical photometry in Fig. [2] 
demonstrates that during the time span between the first and last 
spectroscopic measurements, the continuum flux of GSC 3656-1328 
decreased by a factor of ~ 2, whereas the EWs of Ha, H/3, and H7 
remained constant to within the uncertainties. 
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Fig. 2. — Data and residuals from a microlensing model fit to the 
UV, optical, and IR photometry of GSC 3656-1328. (Top panel) 
The points with error bars show photometry from CBA clear (red) 
and V (cyan), ASAS I (black) and V(blue), Swift M2 (orange), and 
PAIRITEL J (purple), H (grey), and K s (green), as a function of 
JD' =JD-2450000. All data except for the ASAS I and V bands 
have been aligned to the CBA V-band by subtracting a constant 
magnitude offset. For the ASAS data, the additional contribution 
from unresolved neighboring stars has been removed using the mi- 
crolensing model fit. See text for details. The solid line shows 
the best-fit microlensing model. The arrows show the epochs of 
the spectroscopic measurements. The inset shows a detail near the 
peak of the event. (Bottom panel) The residuals from the best-fit 
microlensing model. The achromatic and symmetric variability of 
GSC 3656-1328 is well-fit by a microlensing model. 
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ton noise from the source and comparison stars. In the 
case of the CBA data sets, we estimated the uncertainties 
to be equal to the rms variability about a linear fit to the 
time-series data on the first night ( JD' = 4040) . We align 
all of the bands to the CBA V^-band by subtracting a con- 
stant magnitude offset, except for the AS AS data, which 
due to the larger point-spread function (PSF) and photo- 
metric aperture include an additional contribution from 
unresolved neighboring stars. In this case, we remove 
this unresolved ( "blended" ) light using the microlensing 
model fit, as described below. The baseline magnitude of 
the microlensing source determined in this way is consis- 
tent with that measured directly from higher resolution 
images. As is evident from Figure [2J the variability is 
essentially achromatic (i.e. the color of the event is con- 
stant) over nearly a decade in wavelength. Furthermore, 
although the sampling on the rising side of the event is 
sparse, the variability is apparently symmetric about the 
time of peak brightness near JD' ~ 4040. 

Our spectra show no evidence for any emission fea- 
tures, and are completely consistent with that of a nor- 
mal AOV star. Furthermore, there is no evidence for any 
change in the spectrum of the star between JD' ~ 4043.7 
and 4057.7, during which time the continuum flux de- 
creased by a factor of ~ 2. This can be quantified by the 
EW of the Balmer lines, which were constant to within 
the measurement uncertainties (typically 2% — 4%) dur- 
ing this period (see Fig. []}. Reports of optical spectra 
taken by other groups, which span a broader range of 
epoc hs, apparently conf i rm this lack of spectral e volu- 
tion (|Munari et al.ll2006t iMikolaiewski et al.ll2006allbT ). 

The target-of-opportunity observation by the Swift 
satellite showed no detectable X-ray emission, as would 
be expected from out bursting variables. 

These characteristics of the variability in GSC 3656- 
1328, namely the precise light curve shape including the 
apparent symmetry about the peak, achromaticity, and 
lack of emission features or changes in the absorption fea- 
tures, are expected for microlensing but would be very 
unusual for outbursting variable stars, which generally 
change temperature (and therefore color) during explo- 
sive outbursts. Essentially, what is required to repro- 
duce the type of variability seen in GSC 3656-1328 is 
for the angular size of the star to change by a factor of 
~ y/40 ~ 6 while the temperature remains constant to 
< 5%. The only way for this to occur is for the apparent 
area of the star to change, as occurs in a microlensing 
event. Occultation of the star due to, e.g., an eclipsing 
binary companion can also make such a change in the ap- 
parent area, but obviously in this case one would expect 
a dimming, not the brightening that is observed here. 

When microlensing surveys were first being planned, a 
major concern was the potential contamination from in- 
trinsic stellar variability. As a result, a number of studies 
addressed the question of whether there exists a class of 
variable stars whose variability might be misinterpreted 
as microlensing. At least two types of potential contam- 
inants were identified. 

"Blue bumpers" are blue main-sequence stars whose 
fluxes remain constant for long periods of time, but oc- 
casionally undergo brightenings ("bumps") that are ap- 
proximately symmetric about the peak, and have dura- 
tions that are consistent with the expectations for mi- 



crolensing events (jCook et al.l Il995t lAlcock et al.l 119961 
2000). However, closer inspection reveals that these 
variables have properties that can be used to exclude 
them as the explanation for the variability seen in GSC 
3656-1328. First, many or perhaps all of these vari- 
ables are Be-type stars, exhibiting Balmer line emission 
(|Cook et al.lll995f) . Second, the brightenings are exclu- 
sively low amplitude, with peak brightenings o f less than 
1.5 times the baseline flux (|Alcock et alj |2000). Finally, 
detailed light curves typically show slight asymmetries 
(jCook et alJll995HAlcock et al.ll2000l) . 

At least some dwarf nova (DN) outbursts are charac- 
terized by approximately symmetric brightenings with 
amplitudes and duratio ns that are consistent with that 
seen i n GSC 3656-132 (Idella Valid H99l iBeaulieu et all 
119951: Idella Valle fe Livid 1 19961) . However, it is unlikely 
that the GSC 3656-132 variability is due to a DN out- 
burst. Fi rst, spectra of DNs typically exhibit H or He 
emission (Idella Valid li"99i IBeaulieu et all IT995T) . Sec- 
ond, while the light curves of novae have properties that 
are grossly similar to microlensing e vents, they do not 
follow the standard iPaczvriskl (|1986l ) form at the ~ 2% 
level with whic h the GSC 3656-1328 variability has been 
measured (e.g.. ISmithll2003t lAfonso et al.ll2006fh Finally, 
the interval between outbursts of DNs is known to be 
directly related to the strength of the outburst, such 
that larger outbu rsts generall y have longer interv als be- 
tween outbursts (|Smakl 11984 Ivan Paradiisl |1985[ ). DN 
outbursts with amplitudes similar to the brightening seen 
in GSC 3656-132 would be expected to ha ve average re- 
curre nce times of tens to hundreds of days (|van Paradiisl 
1985]). It seems unlikely that the kinds of observations 
that originally discovered the variability discuss ed here 
(e.g., INakano fe fagol 120061 : iNakano et all 120061 ) would 
have missed previous outbursts; furthermore, there is no 
evidence for such outbursts in the pho tographic plate 
obser vations taken during 1964-1994 (jSamus fc Antipinl 
2006). Regardless, continued monitoring of GSC 3656- 
1328 would certainly allow one to rule out the DN hy- 
pothesis definitively. 

It is interesting to note that the source stars of 
the first two ERO S microlensing event candidates 
(|Aubourg et al.|[l993l ) are both early-type main-sequence 
stars. Spectroscopic observations showed that the source 
of EROS-LMC-1 is a Be star with Balmer emission lines, 
whereas the source of EROS-LMC-2 is a seemingly nor- 
mal AO main-sequence st ar with no apparent emission 
lines (|Beaulieu et al.lll995l ). However, EROS-LMC-2 also 
exhibits periodic variability with an amplitude of ~ 0.5 
mag and a period of ~ 2.8 da ys, indicative of an eclips- 
ing binary (|Ansari et al.lll995h ■ In both cases, continued 
photometric monitoring of the source stars revealed ad- 
ditional brightenings many years later, with amplitudes 
and timescales simi l ar to those of the o riginal events 
(|Lasserre et al.ll2000t[Tisserand et al.ll2007f) . thus exclud- 
ing the microlensing interpretation of the variability. 

Given the evidence in favor of microlensing, we first 
test whether the UV, optical, and IR data can be accept- 
ably fit to a simple microlensing model. This model has 
as parameters the Einstein timescale £e, the impact pa- 
rameter (closest source-lens approach in units of the an- 
gular Einstein radius #e) Uq, and the time-of-maximum 
magnification to. In addition, we fit for the baseline flux 
of the lensed source for each of the eight separate ob- 
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servatory/filter combinations. Since we expect any light 
in the PSF to be completely dominated by the bright 
source, we do not allow for any flux that is blended with 
the source but is not microlensed, with the exception of 
the ASAS data, which are known to contain light from 
unresolved stars in the photometric aperture. Thus, this 
model has 3 + 7 + 2 = 12 parameters. This best-fit model 
is shown in Figure [2] We find that the microlensing 
model provides a reasonable fit to the data: the x 2 /dof 
for the individual data sets ranges from ~ 0.8 for 2802 
data points for the CBA V data set to ~ 2.6 for 34 data 
points in the worst case of the Swift UVM2 dataset. In 
the latter case, the data may be somewhat compromised 
by the fact that the peak of the source PSF was affected 
by coincidence losses in the first two exposures, although 
we attempted to circumvent this difficulty by using an- 
nular apertures. The remaining statistically significant 
deviations from x 2 /dof of unity can be traced to large 
systematic outliers and well-known correlated systematic 
errors in the photometry. In addition, as we discuss in 
detail below, the PAIRITEL data show some evidence for 
a small but non-zero component of blended light which 
may be due to the lens itself. We therefore conclude that 
the microlensing model provides a good representation 
of the data. 

Accepting the microlensing interpretation, we then re- 
fit the microlensing model. We use an iterative procedure 
to remove large outliers from the CBA data (> 3tr for 
CBA V and > 3.5er for CBA clear) and re-normalize the 
uncertainties in each data set by a constant factor to force 
X 2 /dof = 1, except for the UV data, for which we instead 
add a constant 0.03 mag uncertainty in quadrature. To 
provide limits on the magnitude of any contribution of 
light from the lens itself, we now allow for two parameters 
for each of the eight separate observatory/filter combi- 
nations: the flux of the lensed source and a free term 
for any flux that is blended with the source but is not 
microlensed. This model has 3 + 2x8=19 parameters. 
We find t E = 7.19 + 0.03 days, u = 0.0237 + 0.0007, and 
t = 4039.899 ± 0.003. This timescale is unexpectedly 
short; for a typical lens velocity relative to the observer- 
source line of sight of v± ~ 70 km s _1 , a lens mass of 
M ~ 0.3 M Q , and a lens distance of Dl ~ 500 pc, one 
would expect t E ~ 20 days. This implies that the lens is 
moving fast, is of low mass, is very close to the source or 
observer, or some combination of these. 

In order to provide more quantitative constraints on 
the nature of the lens, we perform a Bayesian analy- 
sis. We adopt priors for the parameters of microlensing 
events expected toward the GSC 3656-1328 line of sight 
generated from a Galactic model which includes double 
exponential thin and thick disks, an exponential distribu- 
tion of dust in the vertical direction, and a mass function 
of lenses including stars, brown dwarfs, and remnants. 
We include constraints derived from the measured proper 
motion of the source, as well as constraints from the anal- 
ysis of the light curve. These latter constraints include 
the measurement of the timescale of the event, as well 
as limits on the angular size of the source in units of the 
Einstein ring radius, and the flux of the lens. Details of 
the Bayesian analysis, including information about the 
model assumptions and parameters, as well as the obser- 
vational constraints, are provided in Appendix [Bl 

Figure [3] shows the Bayesian probability densities for 



the properties of the lens star giving rise to the mi- 
crolensing event seen in GSC 3656-1328. Shown are 
the results for the lens mass, distance, relative lens- 
source proper motion fj,, Einstein ring radius Re, and the 
V, J, H, K s magnitudes of the lens assuming it is a main- 
sequence star. The median and 68.3% confidence inter- 
vals are log(M/M ) = — l.Q6±g;l|, D L = 420±l|g pc, 
jti = 43jJ? mas yr^ 1 and R E = 0.35j£^ AU. There is 
an ~ 49% probability that the lens is a main-sequence 
star, an ~ 46% probability that the lens is a brown dwarf, 
and an ~ 5% probability that the lens is a white dwarf. 
The probability that the lens is a neutron star or black 
hole is negligible. The velocity of the lens relative to the 
observer-source line of sight is vj_ = 84Tjs km s _1 , in- 
dicating that it is probably a member of the thick disk. 
Thus the most likely scenario for the lens is that it is a 
low-mass star or brown dwarf in the thick disk, with a 
mass near the hydrogen burning limit. The proper mo- 
tion is likely to be quite high, with (j, > 16 mas yr at 
the 95% confidence level. The apparent magnitudes of 
the lens, assuming it is a main sequence star, are V = 
23.6±l;i, J = 17.3±2-j, H = 16.6±^;|, axidK s = 16.2±^. 
Therefore, light from the lens may be directly detectable 
in a few years, when it has separated from the source, 
using high-resolution, near-infrared imaging. 

Given that the source is quite luminous, and the lens 
likely to be of relatively low mass, one would generally 
expect any blended light due to the microlens to be very 
small by comparison. Indeed the measured blend flux 
values for the CBA I, CBA V, and Swift M2 data sets 
are all consistent with zero and less than 2% of the source 
flux. For the ASAS / and V datasets, the blending is 
significant, but as discussed previously, the ASAS pho- 
tometric aperture is known to contain light from nearby 
stars that are unresolved by ASAS but resolved by the 
CBA data. Surprisingly, we find evidence for significant 
blended light in the J and K s bands, and marginal evi- 
dence in the H band. Specifically, we find blend magni- 
tudes of (J) B = 14.41 ± 0.31, (H) B = 14.76 ± 0.81, and 
(K S ) B = 13.00 ± 0.20. The color and magnitude of this 
blended light are roughly consistent with that of a mid- 
to-late M dwarf located ~ 100 pc away. Although it is 
a priori unlikely that the lens would be sufficiently close 
that its light could be detectable against such a bright 
source, the short timescale of the event already argues 
for a somewhat more nearby lens, making the detection 
of blended light more plausible. Indeed, as can be seen 
from Figure [3l even ignoring any potential constraints 
from the JHK S flux of the lens, there is a small but non- 
negligible probability for the flux from the lens to be as 
large as the measured blended light. 

If we assume that the blended light is indeed due 
to the lens, we can include this information in the 
Bayesian analysis to provide much stronger cons traints 
on the properties of the lens (|Bennett et alj |2006). Fig- 
ure |4] show the resulting probability densities. We find 
that the mass, distance, and proper motion are more 
tightly constrained, log(M/M Q ) = — 0.79±g^|, D L = 

130+4g pc, and jU = 150+20 mas y* -1 ) whereas the con- 
straints on the Einstein ring radius are quite similar 
(Re = 0.39tg;^ AU). The lens apparent magnitudes 

are V = 20.2±i;|, J = U.5t°H, H = l 3 - 9 ±o.i7' and 
K = 13.5j^Q' 2 g. We note that the median posterior value 
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Fig. 3. — Bayesian probability densities for the properties of the 
lens star giving rise to the microlensing event seen in GSC 3656- 
1328. These distributions are derived assuming priors derived from 
models of the mass, velocity, and density distributions of stars in 
the thin and thick disks, and include constraints on the proper 
motion of the source and timescale of the microlensing event, as 
well as limits on the proper motion and V-band magnitude of the 
lens as derived from the light curve. See text. The panels show 
the probabilities for (a) the lens mass (b) the lens distance, (c) 
the lens proper motion, (d) the Einstein ring radius of the lens, 
(e) the V magnitude of the lens, and (f) the J (dotted), H (short 
dashed), and K s (long dashed) magnitude of the lens. In panels 
(a)-(e), the dashed line shows the medians of the distributions, 
and the dark and light shaded regions encompass the 68.3% and 
95.4% confidence intervals, respectively. In panel (e) the dotted 
line shows the 95% c.l. upper limit on the V-band flux from the 
lens. In panel (f), the solid lines show the measured blend fluxes 
from the JHK S light curves. The points show the same along with 
associated uncertainties; the abscissa values are arbitrary. 

of the K s apparent magnitude differs substantially from 
the input constraint due to the strong prior that the lens 
be more distant and less massive, and hence fainter. 

If the blended light is real and indeed due to the lens, 
the lens must be fairly nearby and have a high proper 
motion (/i > 110 mas yr -1 at 95% confidence), and so 
should be resolved from the source in a few years. How- 
ever, there are caveats. The IR data in this event (as 
well as other data taken by PAIRITEL) show evidence 
for low-level systematic errors at the few percent level 
that are correlated on the timescale of several days. Fur- 
thermore, the early IR data were taken when the source 
was sufficiently bright that non-linearity may be impor- 
tant. Either of these systematic effects could easily give 
rise to a spurious blending signal. Regardless, the hy- 
pothesis that the lens is nearby should be testable in the 
near future with high-resolution IR imaging. A measure- 



FlG. 4. — Same as Figure[3] but assuming that the blend fluxes 
measured in the JHK a light curves are real and correspond to flux 
from the lens. This constrains the mass and distance to the lens, 
and implies that the lens is relatively nearby with a large proper 
motion of fj, > 100 mas yr — 1 . Note the greatly compressed scale in 
panel (f) as compared to Figure [3] 

ment of the lens proper motion and relative lens-source 
parallax, when combined with the timescale of the mi- 
crolensing event, would allow for a meas urement of the 
lens mass (|Refsdallll96i lPaczvnski|[l995l) . 

The rate T of microlensing events of a single source ly- 
ing at a distance D by a uniform populat ion of perturbers 
of number density n is (|Paczvhskilll98rl ) 

r = |G 1 / 2 n(m 1 / 2 )M jD 3/2 j (1) 

where G is the gravitational constant, (to 1 / 2 ) is the mean 
square-root mass of the perturbers and (v±) is their 
mean velocity relative to the observer-source line of sight. 
Adopting (to 1 / 2 ) = (0.5 Af ) 1/2 , (v±) = 55kms-\ n = 
0.1 pc -3 , and D = 1 kpc, this yields T = 0.043 Myr" 1 . 
Since there are of order 2 million stars in the Tycho-2 
catalog and noting that the source (GSC 3656-1328) is 
near the brightness limit of this catalog, one would ex- 
pect one microlensing event of a Tycho-2 star every ~ 12 
yr. This crude estimate is in rough agree ment with the 
more detailed calculations of iHanl (|2007f l. who finds a 
rate of one event every ~ 6 yr for stars brighter than 
V = 12. Thus, at first sight, the detection of a mi- 
crolensing event with a source star with V ~ 12 seems 
very plausible, given that amateur and professional ob- 
servers have been combing the skies for comets for almost 
40 years. However, it is unlikely that this microlensing 
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event would have stimulated enough interest to generate 
the high quality follow-up observations that made a con- 
vincing case for microlensing if it had not been magnified 
by at least A > 10. Such events are a factor A^ 1 = 1/10 
less likely to occur. Furthermore, the fraction of the 
event duration when the source is magnified by > A is 
also ~ A^ 1 . Taking these factors into account, we were 
probably lucky to observe such a microlensing event, but 
perhaps not unreasonably so. 

4. DISCUSSION 

The a b ove c alculation and the more detailed study 
by iHanl (|2007f ) show that the event rate for sources 
~ 1 kpc would be quite small, even with a more thor- 
ough and aggressive search procedure that detected es- 
sentially all lenses that come within 1 Einstein radius 
(A ~ 1.34) of the source. However, if the search could be 
extended from 1 to 4 kpc, then the event rate would 
be increased by 4 7 / 2 = 128 times, to roughly 8 per 
year. That is, we obtain D 3 / 2 from equation [T] and 
D 2 from the larger volume probed (since viable targets 
are effectively confined to the two-dimensionsal struc- 
ture of the Galactic plane). Such a survey would en- 
able a probe of the matter distribution near the Sun that 
was equally sensitive to dark and luminous objects (see 
also iDi Stefanol 120051 ) . It would also bring microlens- 
ing planet searches, which have been successful at de- 
tecting some nove l planets against distant sources in 
the Galactic bulge (Bond et al.ll2004 lUdalski et~all l2005l: 
iBeaulieu et alj|2006t iGould et al.ll2006f). nearer to home 
(|Di Stefandl2007ll2008t IDi Stefano fc Night! [200l . 

In order to provide an illustrative example of the planet 
discovery potential of such nearby microlensing events, 
we determine what the sensitivity of this microlensing 
event to planetary companions would have been had it 
been discovered well before peak. Specifically, we create 
a fake data set and then determine the planet detection 
sensitivity of this data set using the methods outlined in 
iDong et alJ (|2006fl . We assume that before JD' = 4034.0, 
the event was sampled at a rate of one point per day. 
For the time intervals of 4034.0 <JD' < 4037.0 and 
4046.0 <JD' < 4080.0, we assume one point every 30 
minutes, and over the peak, 4037.0 <HJD' < 4046.0, we 
assume one point per minute. We assume a photometric 
uncertainty of 0.4% for each data point. These assump- 
tions are reasonable given the large number of amateur 
and professional northern telescopes that are available to 
follow these events, as well as the bright apparent mag- 
nitude of the event. We then determine the planet de- 
tection sensitivity of this simulated data for planet/star 
mass ratios of q= 10" 3 , 10~ 35 , 10" 4 , 10~ 4 ' 5 , and 10~ 5 . 

The finite size of the source will begin to suppress per- 
turbations for mass ratios such that ^fq ~ p„, where 
= 9*/9e is the angular size of the source in units 
of the angular Einstein ring radius 9e- The Bayesian 
analysis presented in JJ] predicts a source size in units of 
the Einstein ring radius of p ~ 10 -2 , and therefore we 
expect perturbations from mass ratios of < 10 -4 to be 
suppressed by finite source effects. For the lowest mass 
ratio we consider, q = 10~ 5 , these suppressions are sig- 
nificant and so we include finite source effects assuming 
p=10- 2 . 

Figure [5] shows the projected positions of the planet 
relative to the position of the primary, where the planet 



would be detected with A% 2 > 160. Here the two com- 
ponents of the position of the planet (b x , b y ) are in units 
of Re, and b x is parallel to the direction of the relative 
source-lens proper motion, such that the source moves 
from left to right relative to the primary along a trajec- 
tory with by = uo = 0.0237. The Bayesian analysis pre- 
dicts a primary mass of M ~ 0.1 M Q and Re ~ 0.4 AU. 
Therefore, a Neptune mass (q ~ 10~ 33 ) planet would 
be detectable over a range of projected separations of 
~ 0.2—1 AU, and planets with masses as low as ~ 0.3 
(q ~ 10 -5 ) would be detectable for some separations near 
the Einstein ring radius (~ 0.4 AU). 
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Fig. 5. — Planetary detection efficiency for simulated data of the 
GSC 3656-1328 microlensing event, assuming that it was monitored 
densely over the peak with a photometric precision of ~ 0.4%. 
The points show the projected coordinates of the planet (b x ,b y ) 
relative to the position of the primary star, at which the the planet 
would be detected with Ax 2 > 160. Here b x is parallel to the 
direction of the relative source-lens proper motion. The projected 
positions arc in units of the Einstein ring radius of the primary; 
the most likely value for the Einstein radius is ~ 0.4 AU. The 
colors correspond to planet/star mass ratios of 10~ 3 (cyan), 10 — 3 5 
(blue), 10 -4 (green), 10 -4,5 (red), and 10 -5 (black). A mass ratio 
of ~ 10~ 4,5 corresponds to an Earth-mass companion for the most 
likely primary mass of ~ 0.1 Mq. 

Identifying all of the microlensing events within 4 kpc 
would require an all-sky (or at least all-Galactic-plane) 
survey that reached a flux level roughly 60 times fainter 
than GSC 3656-1328, or V ~ 16, to compensate partly 
for the greater distance and partly for the increased inter- 
stellar extinction toward more distant targets. It would 
be straightforward to place 120 10 cm lenses, each backed 
by a 20 megapixel camera with 7" pixels, on one single- 
axis mount, and so continuously monitor to the required 
depth the 10,000deg 2 that are within 60° of the zenith 
at any given time. To cover the whole Galactic plane 
would require two such devices, one in each hemisphere. 
Such "fly's eye" telescopes would have heritage in on- 
going experi ments such as ASAS JPojma hski 200^) and 
Super WASP (|Pollacco et al.ll2006D and could be consid- 
ered the "next generation" successors to these exper- 
iments. In addition to detecting nearby microlensing 
events, these telescopes would have a very large num- 
ber of other science applications. Although the technical 
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and survey requirements for each application are fairly 
distinct, these telescopes could in principle be used to de- 
tect thousands of planets as they transit their host star, 
rapidly identify gamma-ray burst afterglows (in plenty 
of time to alert gamma-ray satellites - the reverse of 
the usual procedure), and provide several days' advance 
warning for Tunguska-size meteors that are expected to 
hit the E arth of order onc e per century with 10 megaton 
impacts (lPaczvnskill2006|) . 
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APPENDIX 



PROPERTIES OF THE SOURCE STAR 

GSC 3656-1328 (TYC 3656-1328-1) has a spectral type of A0V-A1V, and is located in the constellation Cassiopeia 
{a = 00 h 09 m 21?99, S = +54°39'43."8 [J2000.0]; I = 116.°8158,6 = -7.°7092). The Tycho catalog (jHflg et al.ll200f)h 
gives a proper motion of fi a = —2.0 ± 2.8 mas yr ,/z<5 = —6.1 ± 2.7 mas yr _1 , and Vr — 11.387 ± 0.076 and 
B T = 11.50 7±0.054. We convert the latter to V = 11.376±0.076 and B = 11.478±0.054 using the standard Hipparcos 
(ESAl ll997f ) transformation. This V magnitude is consistent with our more precise inferred value of V = 11.39 ± 0.01, 
which we adopt here. The color of B — V = 0.102 ± .093 is consistent with, but less accurate than, the value of 
B — V = 0.19 ± 0.01 measured by iMikolaiewski et all ((2006b) when GSC 3656-1328 w as ~ 1 mag above baseline . 
Because we see no evidence for chromaticity, we adopt the more precise determination of iMikolajewski et al.l (j2006b) 
as the baseline color of the source. 

We estimate the distance to the source using its color and magnitude, assuming a dereddenned color and magnitude 
appropriate for its A0V-A1V spectral typ e inferred by IMikolaiewski et al.l (|2006al ). Adopting the spectral type/color 
calibration of iKenvon fc Hartma mi (1995]) we estimate E(B — V)= 0.16 — 0.19, and assuming Ry = 3.1, we estimate 
a U-band extinction of Ay — 0.5 — 0.6. The apparent V magnitude then implies a distance of D = 960 — 1070 pc. For 
definiteness, we adopt Ay = 0.6 and D = 1 kpc. Using the observed V — K s of the source at baseline gives reasonably 
consistent results; however uncertainties about the possible presence of systematics in the PAIRITEL data, as well as 
potential contamination from light due to the lens itself, make these results less secure despite the longer wavelength 
baseline. 

We estimate the angular size of the source to be 6* ~ 10 fxas, ba sed on the der e ddene d (V — K) color and Vq 
magnitude of the source and the color-surface brightness relations of iKervella et al.1 (|2004h . The uncertainty in this 
estimate is not important for our purposes, but it is roughly a few percent due to the uncertainty in (V — K)o and Vq. 

BAYESIAN CONSTRAINTS ON THE LENS STAR PROPERTIES 

In order to provide constraints on the properties of the lens giving rise to the microlcnsing event in GSC 3656-1328, 
we perform a Bayesian analysis, which naturally accounts for priors on the expected populations of microlensing events 
toward the line of sight of GSC 3656-1328, as well as the observed constraints from the light curv e and additional 
(external) information. Our ana l ysis is similar to tha t done for other microlensing events (see, e.g.. lYoo et al J [20041 : 
iDong et al.H2006t iDominikl l2006t Bennett et al.|[2006l ): however, there are some particular nuances in the analysis of 
this particular event that motivate an in-depth discussion. 

We construct prior distributions of microlensing event parameters using simple models of the mass, density, and 
velocity distributions of the lens stars. We adopt double exponen tial models for the thin and thick disks. Our thin 
disk model is the same as that used bv lHan fc Gouldl (|1995l . l2003f ) with a scale length of 3.5 kpc an d scale he i ght o f 
325 pc, and our thick disk model has the same scale length but a scale height of 1 kpc. We adopt a Cha brierl (|2003l ) 



lognormal mass function with a pe ak at M = .079 M w and a width of 0.69 dex for stars with M < M Q and the present 
day mass function as me asured bylReid et al.l (|2002l ) with a logarithmic slope of a = —4.2 for stars with M > M Q . We 
add remnants following Gould (2000;). We adopt Gaussian distributions for the lens velocities that are independent 
of position along the line of sight. These have means in the U,V,W directions of (vjj, Vy, vw) — (0,214,0) km s _1 , 
and one-dimensional velocity dispersions of {ou, cry, aw) — (35, 28, 35) km s _1 for the thin disk and twice this for the 
thick disk. 

In the absence of external constraints, the likelihood of a given lens mass, distance, and velocity, is just the contri- 
bution to the microlcnsing event rate, which for fixed source proper motion and distance is given by, 

£r oc I , T L {p\Gu.\GysGw.\ + PiGu^GyaGw.-i) fx ( -jtj J (Bl) 
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Here subscripts "1" and "2" are quantities for the thin and thick disks, respectively, p is the mass density at the 
position of the lens, and G is the exponential velocity distribution. For example, 



1 ( - \vu -vu 



2 



G " = 7^" P {^1^)- (B2) 

and similarly for the other distributions. The (normalized) mass function of each population of lenses (stars, brown 
dwarfs, etc) is given by (dn/dM)x , and fx is the relative contribution to the number density from each type of lens. 
These are 66.3%, 29.7%, 3.4%, and 0.5% for main-sequence stars+brown dwarfs, white dwarfs, neutron stars, and 
black holes, respectively. The lens velocity relative to the observer-source line of sight is given by 

V_L = Vl,_l - V ,_L I 1 - -p- J - vs,± — . (B3) 

Here vl,±, v S,± and v o,_l are the projected velocities of the lens and the observer perpendicular to observer-source 
line of sight, where vo,_l accounts for both the velocity of the Sun and the velocity of the Earth at the time of the 
event. 

The projected velocity of the source is just its proper motion times its distance, both of which are constrained (§jy. 
We account for these constraints by including a term in the overall likelihood of the form 

£ M oc exp \-\ (»° -J- 2 ' — VI ex P U -J; 6 ' 1 j^Vl ■ (B4) 
M 1 2 V 2.8 mas yr" 1 J j 1 2 \ 2.7 mas yr" 1 J j y ' 

We fix the lens distance at D$ — 1 kpc, but test the effects of changing this value. 

We can account for the constraints from the light curve by including additional likelihood terms. We consider 
constraints from the fitted timescale of the event, the blend flux in the V band, the apparent lack of finite-source 
effects, and, in some cases, the blend fluxes in the JHK S bands. The constraint on the timescale takes the form, 



£t E tx exp 



1 ft E - 7.19 days x L> 



2 V 0.03 days 



(B5) 



The light curve exhibits marginal evidence of finite-source effects, with p* = 0.032 and an upper limit of p* < 0.044 
at the 3cr level. We note that from the observed color and flux of the source, 6* = 10 pas (see $3}, this limit on p* 
constrains the source-lens relative proper motion to be p = #eAe > 10 mas yr" 1 . Since the angular speed of the 
source is 6.4 ± 3.9 mas yr _1 and that of the lens is expected to be similar, this range is plausible. This also implies 
that the lack of finite source effects does not provide a strong constraint on the properties of the lens, other than ruling- 
out lenses that are very close to the source. Nevertheless, we include this constraint on p* from the light curve by first 
determining Axj s (p*), the change in \ 2 from the best-fit point-source model as a function of p*, when allowing all 
the other parameters to vary. The likelihood of a particular p* then takes the form, 

C p , tx exp [-A X ^ s (p*)/2] . (B6) 

We also consider constraints on the flux of the lens arising from measurements of (or limits on) the blend flux from 
the analysis of the light curve. We assume that brown dwarfs and remnants ar e dark. For main -sequence lenses 
we adopt mass-luminosity relations from the solar-metallicity, 1 Gyr isochrone of iSiess et al.l ( 20001). which in turn 
adopt empirical transformations from effective temperature t o standard filter l umin osities from 



Kenvon &: Hartmannl 



(|1995f h We apply a (small) correction to convert the standard iBessell fc Brettl (|1988l) JHK magnitudes given in these 
isochrones to the 2MASS system. 19 To estimate the extinction along the line of sight to the lens, we assume a vertically 
exponential dust disk with a scale height of 120 pc, normalized such that the total extinction to the source is Ay = 0.6 
(see $2}. We use a standard extinction law (Ry = 3.1) to convert to other bandpasses. We then apply the constraints 
on the lens fluxes from the light curve, assuming that the blend flux is entirely due to the lens star. For the y-band 
the strongest constraint on the blend flux is fy = 0.00058 ±0.00181, where the units are such that fy = 1 corresponds 
to a V — 10 star. If we are considering only the constraint from the F-band flux, then the likelihood is simply 

C lc ex exp (-Axl/2) , (B7) 

where A\ 2 = [{fv — fv, model )/°v] 2 ; °V = 0.00181, and jV, mo dci is the V-band flux of the lens predicted by the 
model. We can also include constraints from the blend fluxes in the infrared. The measured blend fluxes fj,fn,fK B 
in the JHK S bands are given in $3] Since these fluxes are correlated with each other and with the y-band flux 
within the microlensing fit, we must adopt a somewhat more sophisticated approach. We first construct the vector 
= (fv - ,/V,modci,/j - /j,model)/ff - /ff,model) fj( a - fx s , model) , where /y,model, /j.model, ••• are the model fluxes. 
From the microlensing fit to the light curve, we can also construct the covariance matrix C for the four blend-flux 
parameters. The difference in x 2 between the fluxes predicted by the model and the measured fluxes is 

4 4 

«=1 j=\ 

19 See |http://www.astro. caltech.edu/~jmc/ 2mass/v3/transformations/ 
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where B = C 1 . The likelihood is then determined using Equation (|B7|) , as before. 
Finally, the likelihood of a particular parameter combination is given by 

£*ot = £r£ti,£tE£p t £ic- (B9) 

We construct a posteriori probability densities using the Markov chain Monte Carlo method. We first randomly 
choose values for the lens mass, distance, and (U, V, W) velocity components vu,vv,v\y, over a range of parameter 
space that is broad in comparison to the posterior probability distributions. We also randomly choose whether the 
lens is a main-sequence star, brown dwarf, white dwarf, neutron star, or black hole. Finally, we choose a random 
value for the two components of the proper motion of the source. We evaluate the relative likelihood of this parameter 
combination using equation IB9I We then randomly move to another point in the parameter space of lens mass, 
distance, velocity, and source proper motion, and evaluate the likelihood of this new parameter combination. We step 
in parameter space by adding to each of the parameters a random value drawn from a Gaussian distribution with zero 
mean, and dispersion chosen to efficiently sample the likelihood surface. Specifically, these dispersions are 0.043 dex 
in log(M/M Q ) for M, 148 pc for Dl, ajj,crv,o-w for the lens velocities, and 2.7 mas yr _1 for the two components 
of the source proper motion. If the ratio of the likelihood of the new parameter combination to the old combination 
is greater than unity, we take the step. Otherwise, we draw a random number between and 1. If this number is 
less than this ratio of likelihoods, we take the step; otherwise it is rejected. After discarding the first 10 4 steps, we 
record every 10 3 steps until w e form a chain of 10 5 poi nts. We form 10 such chains, each starting from different initial 
conditions, and calculate the iGelman fc Rub in (1992) i?-statistic. This is within 0.2% of unity for each parameter, 
indicating that the chains are well-mixed and converged. We then merge the chains and use the result for the final 
probability distributions. 

We also test the effects of changing the source distance by ±100 pc and the F-band extinction to the source by 
±0.2 mag. We find that the medians of the probability distributions change by < 5% for all of the parameters of 
interest. 
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